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Transport of molecules across membrane channels is investigated theoretically us- 
ing exactly solvable one-dimensional discrete-state stochastic models. An interaction 
between molecules and membrane pores is modeled via a set of binding sites with 
different energies. It is shown that the interaction potential strongly influences the 
particle currents as well as fluctuations in the number of translocated molecules. 
For small concentration gradients the attractive sites lead to largest currents and 
fluctuations, while the repulsive interactions yield the largest fluxes and dispersions 
for large concentration gradients. Interaction energies that lead to maximal cur- 
rents and maximal fluctuations are the same only for locally symmetric potentials, 
while they differ for the locally asymmetric potentials. The conditions for the most 
optimal translocation transport with maximal current and minimal dispersion are 
discussed. It is argued that in this case the interaction strength is independent of 
local symmetry of the potential of mean forces. In addition, the effect of the global 
asymmetry of the interaction potential is investigated, and it is shown that it also 
strongly affects the particle translocation dynamics. These phenomena can be ex- 
plained by analyzing the details of the particle entering and leaving the binding sites 
in the channel. 
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I. INTRODUCTION 

Membrane channels are large water-filled hollow protein structures that control the 
transport of metabolite molecules between different cells or between different cellular 
compartments.-"^ These processes are critically important for biological systems, and re- 
cent experimental and computational studies suggest that, contrary to earlier views, perme- 
ation of molecules across such large channels is efficient and selective.-"^'^'^'^'^'^ However, our 
understanding of these phenomena is still very limited. 

To analyze mechanisms of efficiency and selectivity of the transport across membrane 
channels several theoretical methods have been presented.— >^>^>^>^>^>^ One approach uti- 
lizes a continuum description where a single-molecule transport across membrane pores is 
viewed as an effective one-dimensional diffusion in a potential of mean forces created by 
interactions between the solute and protein channel.— '^^'^'i^ The interactions are modeled 
as square well potentials that occupy the entire channel. Using this method, it was shown 
that the particle's current can be maximized for some interaction strength that depends on 
the solute concentrations, diffusion constants and geometry of the pore.— It is also possible 
to compute the inter-channel potential of mean force that maximizes the flux.— We recently 
developed a discrete-state stochastic model of the channel-facilitated membrane transport 
that takes into account the existence of binding sites inside the pore.— By mapping the 
discrete-state model of the permeation through the pore to a single-particle hopping along 
a periodic lattice, the particle currents have been obtained explicitly for all sets of param- 
eters. The theoretical analysis suggested that the presence of the binding sites accelerates 
the particle flux for small concentration gradients, while the repulsive binding sites are more 
advantageous for creating the most optimal current for large concentration differences. In 
addition, it was shown that the asymmetry in the interaction potential, e.g., the spatial 
position of the binding sites, might also significantly change the particle dynamics. Similar 
observations have also been obtained in the continuum models of membrane transport.—*^ 
Theoretical calculations show that both continuum and discrete-state descriptions are closely 
related,—*^ and the results obtained by these approaches can be mapped into each other. 

In many biological systems concentrations of molecules that involved in permeation 
trough the membrane pores are rather small, and this points out to the importance of 
fluctuations in the number of translocated particles.^ However, current theoretical studies 
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concentrate mostly on the description of fluxes, i.e., the average number of particles moved 
across a single pore per unit time.—^^^^^^^'^^^^^ In this work, we analyze the effect of in- 
teractions on fluctuations in the number of translocating molecules and compare it with 
the effect of interaction potentials on the particle currents. Our approach is based on the 
discrete-state stochastic models for which all dynamic properties can be calculated explicitly 
via mapping to the single-particle random hopping model on periodic lattices.— 



II. MODEL 



We consider a transport of particles through a membrane channel as an effective one- 
dimensional motion across a cylindrical pore with binding sites inside. The pore separates 
two chambers as illustrated in Fig. 1. The molecule's concentrations in the left and right 
chambers are ci and C2, respectively. The concentration gradient Ac = ci — C2 (ci > C2) 
drives the particle current mostly from the left to the right. The molecule can enter the 
channel from the left with the rate uq = kon.ci, but from the flrst binding site it can also 
return back with the rate Wi = k^ff. Similarly, the molecule can jump to the pore from the 
right chamber or jump out from the state with rates wat+i = Wq = fconC2 and u^r = ^o//j 
correspondingly. The particle at the site j {j = 1,2, ■■■ , A^) can move one site forward 
with the rate uj, or it can jump backward one site with the rate Wj: see Fig. 1. Because 
the concentration of solute molecules is typically small, it is assumed that particles do not 
interact with each other, and there is no more than one particle can be found in the channel 
at all times. The probability to flnd the particle at the binding site j at time t is given by 
a function Pj(t), and time evolution of the translocation process can be described by set of 
master equations, 

dm 



u 



i-lPj-l{t) + Wj+iP,+i{t) - {u, + W,)Pj{t), (1) 



dt 

where j = 1, 2, ■ ■ ■ , A^, and we have deflned 

N 

PoW = P7v+i(t) = l-l]P,(t) (2) 

i=i 

as the probability of flnding the channel empty and the molecule outside of the pore at time 
t. 
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Each binding site corresponds to a minimum in the free energy profile for translocation 
as illustrated in Fig. 2. We associate the strength of the interaction at site j with the 
parameter {—£j), assuming that zero free energy is at the entrance, Eq = 0.— Note that 
the free energy at the exit to the right chamber is {—En+i)- Then, Ej > describe the 
attractive binding sites (with respect to the left chamber), while negative Ej corresponds to 
the repulsive sites.— The transition rates between the sites are related to binding energies 
via detailed balance conditions, 

Uj{Ej,Ej+i) _ Uj{Ej,Ej+i = 0) _ Uj{Ej = 0,Ej+i'' 



,X,+i = -^ ' (l/x^, (3) 

Wj+i[Ej,Ej+i) Wj+i(£j,£j+i = 0) Wj+i[Ej = Q,Ej+i) 

with Xj = exp [Ej/kBT]. Dynamic properties of the system depend on explicit expressions 
for the transition rates that can be written in the following form,— 

Uj{Ej+i) = Uj{Ej+i = 0)x/, Wj+i{Ej+i) = Wj+i{Ej+i = 0)x/ , (4) 

or 

= Uj{Ej = 0)X^'~\ Wj+i{Ej) = Wj+i{Ej = 0)x^/, (5) 

where < 6j < 1 are interaction-distribution coefficients that describe how binding en- 
ergies are distributed between forward and backward transition states. These coefficients 
also provide relative distances between neighboring free-energy minima and transition states 
(maxima in Fig. 2). Similar parameters have been utilized in the analysis of motor pro- 
tein's dynamics. To simplify calculations, in this paper we will assume that interaction- 
distribution coefficients are the same for all binding sites, i.e., 6j = 6 for all j. 

If we define M{t) as number of particles that translocated through the membrane channel 
at time t, then the stationary-state particle flux is given by 

(6) 

where averaging is taken over all possible translocation events. To specify fluctuations, a 
dispersion D is introduced in the following way, 

2t-oo dt ^ ' 

The discrete-state stochastic model of channel-facilitated membrane transport can be 
solved exactly at large times by utilizing the mapping to the single-particle hopping model 
along infinite one-dimensional periodic chain.— This mapping can be understood using the 
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following arguments. Consider multiple identical membrane channels (as shown in Fig. 1) 
arranged in a sequence such that the particle exited from one channel can enter the next 
one. At the stationary state the flux is constant, and the transport of the particles along 
the sequence of channels (with binding sites in each) is identical to the motion of the 
single particle along one-dimensional periodic lattice with a period of + 1 sites. The 
number of states in the period of the effective lattice is larger than in the channel because 
the additional state corresponds to the situation when the particle is outside of the channel. 
For the effective single-particle hopping model all dynamic properties are known exactly.— 
Thus, transport across membrane pores can be analyzed explicitly for all sets of parameter 
in the stationary-state limit. 



III. RESULTS AND DISCUSSION 



To investigate the effect of interactions on particle dynamics in channel-facilitated mem- 
brane transport we consider the simplest model with A^ = 1 binding site inside the pore.— 
The binding energy is equal to —e, and the detailed balance conditions can be written as 

uo{e) UQ{e = 0) ui{e) ui{e = 0) 



-X, 



wi{e) wi{e = 0) ' wo{e) wo{e = 0) 
where x = exp [e/ksT]. The transition rates are given by^ 



(8) 



Uo{e) = uqx^, Wi{e) = Wix^ ^, Ui{e) = UqX^ ^, Wo{e) = WqX^. 
Then using known results,— the particle flux for A^ = 1 model has a simple form, 

{UqUi - WoWi)x^ kon{Ci - C2)x^ 



(9) 



J 



{uq + Wq)x + {ui + Wi) 2 + ^22i(£l±£2) 



(10) 



•'off 



X 



The corresponding expression for the dispersion is more complex 
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D = - 

2 [uq + Wojx + ui+wi 2 

It is convenient to consider the relative flux, 



kon{Cl + C2)X + 2koff 



J 

To 



koniCl + C2) + 2ko 



ff 



kon{Cl + C2)X + 2koff 



(11) 



(12) 
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and the relative dispersion, 



D 
Do 



koJci + C2) + 2k. 



off 



konici + C2)x + 2k, 



off 



h h 



onKoff[Cl + C2) 



X 



, fconfcg//(ci-C2)^ 
' [konici+C2)x + 2kof f] 



-X 



konkoffiCl + C2) 



fcLfco//(ci-C2)^ 



(13) 



where Jo and -Do are dynamic properties of the system without interactions {e = 0). 

The molecular flux and dispersion are strongly influenced by interactions at the binding 
site, as illustrated in Fig. 3. The presence of strongly attractive or strongly repulsive 
binding sites lead to decrease in both particle currents and fluctuations. For intermediate 
interactions the molecular flux and dispersion are large, and this behavior is independent 
of concentration gradients. The relative particle current reaches a maximum value at the 
interaction energy Ej that can be obtained from Eq. (|T2l) . 

e 2koff 



ksT In 



(14) 



1-9 kon{Cl + C2) 

Fluctuations in the number of translocating particle also produce a maximum as a function 
of the binding energy strength, as can be seen in Fig. 3, but the corresponding energy of 
interactions e]-, yields a more complex expression, namely, 

koff 



G 



. 1 - 6* / konCl _ 

for the simplest case of C2 = 0. The auxiliary function G{'y) is deflned as 



(15) 



3 64 , 4 



51274 - 97672 + 512 4 

1 — 7. 

27 3 ' 



(16) 



It has the following properties that G(7 = 1) = 2, and for small 7 we have ^(7) ^ 7/3, 
while for 7^1 the asymptotic behavior is different, ^(7) ~ 47. 

Interactions energies e} and that produce maximal fluxes and maximal dispersions 
depend on the concentrations outside of the membrane channel and on the local environ- 
ment around the binding site (via the interaction-distribution parameter 6) as shown in Fig. 
4. When the concentration of particles in the left chamber is small the presence of attrac- 
tive site leads to the largest particle current, and at these conditions fluctuations are also 
maximal. However, for large concentrations Ci the repulsive site produces the biggest flux 
and dispersion. This behavior has been explained before by considering the details of the 
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particle dynamics near the binding site.— It can be shown that the total time to move across 
the channel consists of two terms that describe the effective time to enter and to leave the 
binding site.— The maximal current is achieved when these two terms are of the same order. 
Then for small concentration gradients the effective time to enter the binding site is large, 
and to produce the optimal flux it is required that the particle stays longer in the channel, 
which corresponds to > 0. For large values of C\ the time to move into the binding site 
is small, and only repulsive interactions will lower the time for the particle in the channel, 
producing the largest current. It is reasonable to suggest that similar arguments can be 
used to explain the behavior of dispersion. 

The results presented in Fig. 4b also show that interaction energies that lead to maximal 
fluxes and dispersions strongly depend on the location of transition states near the binding 
site. When the position of the transition state between two minima in the potential of mean 
forces is closer to the right one [Q > 0.5), the binding interactions that produce maximal 
fluctuations are larger than the interactions that yield the maximal flux. At the same time, 
for the transition state closer to the left state {6 < 0.5) the situation is different, and Ej > e*j^. 
Thus for the locally symmetric potential of mean forces {6 = 0.5) the interaction energies 
e*j and e}j coincide, while the local asymmetry {6 ^ 0.5) yields different values for the 
most optimal binding sites interactions. This observation can be understood by analyzing 
Eqs.f ll2p and f|T3|) . Comparing these two expressions for the maximal interactions strengths 
it can be shown that 

where the function F{x), the ratio of the relative dispersion over the relative flux, is given 
by 

KonKoff[Ci + C2) / [fc„„(ci+C2)+2fc„^^P 

Taking derivative of the left and right side of Eq. (fT7|) with respect to the variable x, we 
obtain 

d{D/Do) ^ 
dx 

~xkon{Cl+C2) 



dx 
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Then '^(-^/^o) and '^^'Y'^"^ simultaneously become equal to zero only if 

ax ax ^ T. J 

x=-p^,. (20) 

kon[Cl + C2) 

Comparing this result with Eg. (1141) and recalling that x = exp {e/ksT) it can be concluded 
that e*j = e}j only for 0=1/2. 

The observation that the maximal particle current and the maximal dispersion for 6 7^ 0.5 
are realized for different interaction strengths raises the question of what is the binding 
interaction energy that allows us to have the largest possible current simultaneously with 
smallest fluctuations in the number of translocated molecules. From the point of view of 
functionality of cellular processes this interaction strength might be viewed as the most 
optimal. To answer this question we analyze the function F{x), the ratio of the relative 
dispersion over the relative current, as given in Eq. (ITSl) . The most optimal conditions can 
be reached by minimizing this function. It can be shown that this ratio is minimal when 

2koff 



^Ipt = ksT In 



konici + C2) 



(21) 



Note that, as discussed above, for the locally symmetric potentials {6 = 0.5) we have = 
e}) = elp^. However, generally the most optimal interaction does not produce the largest 
fluxes or fluctuations, and it is independent of the interaction-distribution parameter 6. 
These results are illustrated in Fig. 5. 

Theoretical and experimental investigations of the potential of mean forces^"^*^ indicate 
that the free-energy landscape for the molecules permeating across membrane channels is 
generally globally asymmetric. This asymmetry is important for biological channels,— 
and for the transport across artificial pores. It was shown recently by one of us-^ that the 
asymmetry influences the particle fluxes, and the origin of this phenomenon was discussed 
by analyzing the dynamics of particle translocation. To study the effect of asymmetry on 
fluctuations in the number of translocating particle, we consider two membrane channel 
models with N = 2 binding sites. In the first model, the binding energies on sites 1 and 2 
are equal to {—e) and 0, correspondingly, while in the second model the order is reversed. 
Otherwise both models are identical. Putting the special binding site [with energy {—e)] on 
the first or on the second site introduces the asymmetry in the system. For the channel in 
the first model the particle current is equal to^i^^ 

koniCl - C2) 



Jl 



konCl _|_ -0 I 1 _|_ fcp// _|_ fconC2 _i_ fconC2 \ _i_ ryl-d ( fconC2 _l_ fcpnCl \ _i_ „konC± 



(22) 



9 



where we assumed that transition rates inside the channels are the same, Ui = W2 = ol. For 
the model with the interaction in the second binding site one finds 

T fco«(ci - C2) , . 

J2 — T 7 7 N 7 N T- 1^<JJ 

''off V " ''°ff " / V ''"ff " / . 

The explicit expressions for dispersions Di and D2 for both models can be written,— 
however they are quite bulky and we will not present them here. To measure the effect of 
the asymmetry on dynamic properties we plot the ratio of currents, J1/J2, and the ratio 
of dispersions, D1/D2, for both models in Fig. 6. For all quantities deviations from unity 
indicate that the asymmetry is important for particle currents and for fluctuations in the 
number of translocated molecules. For all concentrations outside of the membrane channel 
the repulsive interaction {e < 0) on the flrst binding site leads to larger particle current and 
dispersion, although fluctuations are affected less strongly, in comparison with the situation 
when the repulsive interaction is on the second binding site. However, putting the attractive 
site flrst {e > 0) have an opposite effect: the particle currents and dispersions are generally 
lower for the flrst model. 

To understand these phenomena let us consider channels with repulsive interactions. In 
the flrst model, with the repulsion on the flrst binding site, the particle spends more time in 
the second site because this position is energetically more favorable for the particle. Then 
the molecule is closer to the right chamber, and it might easily exit out. This leads to 
larger fluxes and dispersions. In the second model, with the repulsive interaction on the 
second site, the particle stays longer in the flrst binding site, which is further from the exit, 
producing smaller currents and fluctuations. These arguments suggest that in the transport 
of molecules across the membrane channels it is possible to control molecular fluxes without 
influencing much the fluctuations by putting the binding sites at the proper positions. It also 
agrees with the idea of the optimal inter-channel potential developed in continuum models 
of channel-facilitated transport.— 

IV. SUMMARY AND CONCLUSIONS 

We presented a theoretical investigation of the transport of molecules through membrane 
pores by analyzing discrete-state stochastic models that allow one to calculate explicitly 
dynamic properties of the system. It was shown that interaction potentials between the 
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molecule and the channel strongly affect translocation dynamics. For small concentrations 
outside of the membrane the attractive binding sites produce largest particle currents and 
dispersions, while for large concentrations the repulsion leads to large fluxes and fluctua- 
tions. For locally asymmetric potentials (with the interaction-distribution coefficient 9 not 
equal to 0.5) maxima in the particle currents and dispersions are achieved for different in- 
teraction strengths, and for locally symmetric potentials {0 — 0.5) the largest molecular 
fluxes and fluctuations are taking place for the same interactions. We found conditions for 
the most optimal transport across the membrane pores, when the largest possible flux is 
accompanied by the smallest possible fluctuations. It was shown that the most optimal 
interaction strength is independent of the local asymmetry in the potential, although it gen- 
erally does not coincide with the positions of maximal fluxes or maximal dispersions. We 
also investigated the effect of the global asymmetry on translocation dynamics, and it was 
argued that the location of binding sites with different interaction strengths strongly affects 
the molecular transport across the channels. These phenomena are explained by using the 
details of dynamics of particle entering and leaving the binding sites. Our theoretical analy- 
sis suggests a possible mechanism of selectivity and efficiency of membrane channels: tuning 
the interaction potential by changing the interactions and asymmetry in the potential (both 
local and global) it is possible to control the translocation dynamics. 

It is important to note that our theoretical approach is based on several oversimplifled 
assumptions. Speciflcally, interactions between the molecules and three-dimensional nature 
of membrane channels and interaction potentials are neglected. It will be important to com- 
pare our theoretical predictions with more realistic theoretical models and with extended 
molecular dynamics computer simulations. However, the most important test of our theoret- 
ical approach should come from the experiments that will simultaneously measure molecular 
fluxes and fluctuations. We believe that a combination of analytical, computation and ex- 
perimental methods will help to uncover the mechanisms of translocation across membrane 
channels. 
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Figure Captions: 

Fig. 1. General schematic view for discrete-state stochastic models of channel-facilitated 
transport. A cylindrical membrane divides the system into three parts: the left chamber 
with particle concentration Ci, the right chamber with particle concentration C2, and the 
pore which can be occupied by a single particle. Open circles correspond to the binding 
sites in the channel. At site j the particle jump forward and backward with rates Uj and 
Wj, respectively. The filled circle describes the position currently occupied by the particle. 

Fig. 2 Potential of mean forces for the channel-facilitated membrane transport. The free 
energy at the entrance (site 0) is equal to zero. Sites 1 and 2 are attractive, while the site 3 
is repulsive. The site + 1 corresponds to the right chamber. 

Fig. 3. Relative molecular fluxes and relative dispersions as a function of the interaction 
strength for the model with A^ = 1 binding site for different concentrations and for different 
interaction-distribution factors. The transitions rates, kon = 15 fiM^^s^^ and koff = 500 
s^^, are taken from Ref.- For all calculations C2 = is assumed, a) ci = 10 /iM; b) ci = 500 
/iM. 

Fig. 4. Interactions producing maximal current and dispersions as a function of a) the 
external molecular concentration ci, and b) the interaction-distribution parameter 6. The 
transitions rates, kon = 15 ;uM~^s^^ and koff = 500 s~^, are taken from Ref.- For all 
calculations C2 = is assumed. 

Fig. 5. The ratio of relative dispersion over the relative current as a function of the in- 
teraction strength for the model with A^ = 1 binding site for different concentrations. The 
transitions rates, kon = 15 fiM~^s~^ and koff = 500 s^^, are taken from Ref.- For all calcu- 
lations C2 = is assumed. 

Fig. 6. The ratio of current and dispersions as a function of the interaction strength for two 
models with N = 2 binding sites. The transitions rates, kon = 15 fiM^^s"^ and koff = 500 
s~^, are taken from Ref.~ For all calculations C2 = 0, 6* = 0.5 and a = koff are assumed. 
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Figure 1. Kolomeisky and Kotsev 
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Figure 2. Kolomeisky and Kotsev 
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Figure 3a. Kolomeisky and Kotsev 
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Figure 3b. Kolomeisky and Kotsev 
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Figure 4a. Kolomeisky and Kotsev 
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Figure 5. Kolomeisky and Kotsev 




Figure 6. Kolomeisky and Kotsev 



